source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days truncated.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days truncated.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days truncated.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days truncated.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days truncated.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days truncated.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days truncated.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
View(devel_time)
View(devel_time)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
# Descriptive statistics
attach(devel_time)
describe(male)
describe(female)
detach(devel_time)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
p*(1-p)^(t-1)
t*p
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
p <- 0.5
t <- 1:9
y <- sum(t*p*(1-p)^(t-1))
y
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
p <- 0.15
t <- 1:9
y <- sum(t*p*(1-p)^(t-1))+10*(1-p)^9
y
p + 2*p*(1-p) + 3*p*(1-p)^2 + 4*(1-p)^3 + 5*p*(1-p)^4 + 6*p*(1-p)^5 +
7*p*(1-p)^6 + 8*p*(1-p)^7 + 9*p*(1-p)^8 + 10*(1-p)^9
t*p*(1-p)^(t-1)
9*p*(1-p)^8
8*p*(1-p)^7
7*p*(1-p)^6
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
t*p*(1-p)^(t-1)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity ptoids.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity ptoids.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
sum(p*(1-p)^(t-1))
(1-p)^(maxt-1)
sum(p*(1-p)^(t-1))+(1-p)^(maxt-1)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity ptoids.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
View(outdat12)
View(outdat12)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
View(outdat12)
View(outdat12)
outdat12
outdat12[84,]
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity ptoids.R', echo=TRUE)
tables(devel_time)
table(devel_time)
?table
table(females)
attach(devel_time)
tables(females)
detach(devel_time)
attach(devel_time)
table(females)
detach(devel_time)
attach(devel_time)
table(female)
detach(devel_time)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity weevils.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Findp.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Control simulation model days geometric.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity ptoids.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/Adult longevity ptoids.R', echo=TRUE)
?prod
nsurv <- n0*cumprod(1-mu)
mu <- c(0,0,0,0,0,0,0.15,0.25,0.50,0.75,1)
n0 <- 100
nsurv <- n0*cumprod(1-mu)
nsurv
nage <- nsurv[2:10]-nsurv[1:10]
nage <- nsurv[2:10]-nsurv[1:9]
nage
nage <- nsurv[1:9]-nsurv[2:10]
nage
nsurv
nage
nsurv
length(nsurv)
nsurv <- n0*cumprod(1-mu)
nage <- nsurv[1:10]-nsurv[2:11]
nage
ybar <- sum((1:10)*nage)
ybar
ybar <- sum((1:10)*nage)/n0
ybar
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
mu
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
source('C:/Users/John/Desktop/Weevils/Dishes-days/mu_schedule.R', echo=TRUE)
?apply
?sapply
source('~/My Documents/Biostats 17/Chapter 7 Code/R Code/Linear.R', echo=TRUE)
?rpois
?reshape
?matrix
y5 <- rpois(5*n,n,5)
y5 <- rpois(5*n,5)
y5 <- matrix(y5,nrow)
y5 <- matrix(y5,n,5)
y5 <- matrix(y5,n,5)
View(y5)
View(y5)
y5 <- apply(y5,1,mean)
source('~/My Documents/Biostats 17/Chapter 7 Code/R Code/central_limit_theorem.R', echo=TRUE)
?sapply
?apply
source('~/My Documents/Biostats 17/Chapter 7 Code/R Code/Linear.R', echo=TRUE)
?sapply
source('~/My Documents/Biostats 17/Chapter 8 Code/R Code/Likepois_random_np.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 8 Code/R Code/Likepois_random_np.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 8 Code/R Code/Likepois_random.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 8 Code/R Code/Likenorm_random.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 9 Code/R Code/confidence_intervals.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 10 Code/R Code/one-sample_t_test_power2.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 11 Code/R Code/bark_beetle_experiment.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 11 Code/R Code/bark_beetle_random.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 11 Code/R Code/bark_beetle_experiment_ttest.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 11 Code/R Code/bark_beetle_experiment_ttest.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 12 Code/R Code/oneway_power.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 12 Code/R Code/oneway_power2.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 12 Code/R Code/oneway_power3.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 13 Code/R Code/Kneitel_2010_algae_lsd2.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 13 Code/R Code/Kneitel_2010_algae_tukey2.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 13 Code/R Code/Kneitel_2010_algae_tukey2.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 13 Code/R Code/Kneitel_2010_cover2.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 13 Code/R Code/Kneitel_2010_algae_fdr2.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Maestre_biomass_new.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Maestre_biomass_new.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Tdubius_bait_color_new.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
?lsmeansLT
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Maestre_biomass_new_unbalanced.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Maestre_biomass_new_unbalanced.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
install.packages("RLRsim")
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
# TrapRCBD_clerids.R - trapping experiment data
# Randomized block mixed model and Tukey method
# Load necessary libraries
library(ggplot2)
library(lme4)
library(lmerTest)
library(multcomp)
# Read in data set
trapdata <- read.table(header=T,colClasses=c("factor","factor","numeric"),text="
block treat count
1	AP	    4
1	BLANK   0
1	FRAP	 79
1	IDAP	  7
1	ISAP	 10
2	AP	    1
2	BLANK	  0
2	FRAP  124
2	IDAP	 13
2	ISAP	 20
3	AP	    0
3	BLANK	  0
3	FRAP	 14
3	IDAP	 NA
3	ISAP	  2
4	AP	    0
4	BLANK	  0
4	FRAP   15
4	IDAP	 11
4	ISAP	  7
5	AP	    0
5	BLANK	  0
5	FRAP	 29
5	IDAP	  7
5	ISAP	  7
6	AP	    2
6	BLANK	  0
6	FRAP	 70
6	IDAP	 14
6	ISAP	 20
")
# Apply transformations here
trapdata <- transform(trapdata,y=log(count+1))
# Print data
trapdata
# Delete blank traps before analysis
trapdata2 <- subset(trapdata,subset=treat!="BLANK")
# Plot treatment means and confidence intervals
ggplot(trapdata,aes(treat,y))+
geom_jitter(size=2,position=position_jitter(width=0.1))+
stat_summary(fun.y="mean",geom="point",size=3)+
stat_summary(fun.data="mean_cl_normal",geom="errorbar",width=0.2,size=1)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
?exactLRT
library(nlme);
data(Orthodont);
##test for Sex:Age interaction and Subject-Intercept
mA<-lme(distance ~ Sex * I(age - 11), random = ~ 1| Subject,
data = Orthodont, method = "ML")
m0<-lm(distance ~ Sex + I(age - 11), data = Orthodont)
summary(mA)
summary(m0)
exactLRT(m = mA, m0 = m0)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
View(lmH0)
View(lmH0)
?lmer
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
anova(lmeH1,lmH0)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
?exactLRT
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
?lmer
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
AIC(lmH0)
ls(lmH0)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
?anova
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
drop(lmeH1,(1|block))
drop(lmeH1,block)
drop(lmeH1,scope=(1|block))
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
?lm
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/TrapRCBD_clerids_block_test.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Rand_treatments.R', echo=TRUE)
?arrange
?order
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Rand_treatments.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Rand_treatments.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Rand_treatments.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Rand_treatments.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 14 Code/R Code/Rand_treatments.R', echo=TRUE)
summarize(tttdata)
source('~/My Documents/Biostats 17/Chapter 15 Code/R Code/arcsine.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 15 Code/R Code/TrapRCBD_clerids.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 16 Code/R Code/WKWtest_chitons_Venados.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 16 Code/R Code/WKWtest_chitons_Venados.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 16 Code/R Code/WKWtest_chitons_Venados.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 16 Code/R Code/KWtest_chitons_3islands.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 16 Code/R Code/KStest_chitons_Venados.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 16 Code/R Code/Randtest_bobcat_parasites.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 17 Code/R Code/SPBattack.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 17 Code/R Code/logistic.R', echo=TRUE)
anova(regout)
source('~/My Documents/Biostats 17/Chapter 17 Code/R Code/SPBattack.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 17 Code/R Code/logistic.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 17 Code/R Code/SAprob2.R', echo=TRUE)
?cor.test
?cor.prob
setwd("~/My Documents/Biostats 17/Chapter 18 Code/R Code")
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/Iris.R', echo=TRUE)
cor.test(irisdata)
library(Hmisc)
rcorr(irisdata)
rcorr(as.matrix(irisdata))
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/Iris.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/Iris.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/Iris_all.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/rhocalc.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/Iris.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/rhocalc.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 18 Code/R Code/rhocalc.R', echo=TRUE)
setwd("~/My Documents/Biostats 17/Chapter 19 Code/R Code")
setwd("~/My Documents/Biostats 17/Chapter 19 Code/R Code")
source('~/My Documents/Biostats 17/Chapter 19 Code/R Code/Nested_ANOVA_Anagrus.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 19 Code/R Code/Maestre_biomass_3way_new.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 19 Code/R Code/ANCOVA_fitness.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 20 Code/R Code/gof_peas.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 20 Code/R Code/gof_peas.R', echo=TRUE)
install.packages("DescTools")
source('~/My Documents/Biostats 17/Chapter 20 Code/R Code/gof_peas.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 20 Code/R Code/age_structure.R', echo=TRUE)
CrossTable
?CrossTable
source('~/My Documents/Biostats 17/Chapter 20 Code/R Code/age_structure.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 20 Code/R Code/age_structure.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 20 Code/R Code/chytrid.R', echo=TRUE)
source('~/My Documents/Biostats 17/Chapter 19 Code/R Code/Nested_ANOVA_Anagrus.R', echo=TRUE)
